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1.  Introduction 

There  has  been  extensive  study,  both  experimentally  and  theoretically,  on  how  a 
semiconductor  nanocrystal  evolves  to  a  bulk  material  and  how  its  corresponding  electronic 
properties  change.  *  Experimental  observations  show  a  close  relationship  between  the  size  of 
nanocrystals,  the  lowest  energy  optical  absorption  peak  and  the  highest  energy  band  of  the 
luminescence  spectrum.^  These  are  usually  assigned  to  the  transition  between  the  valence  band 
and  the  conduction  band  or,  in  the  molecular  analogy,  between  HOMO  and  LUMO.  Among  the 
significant  issues  are  the  size-dependent  shift  of  the  lowest  energy  transition  and  its  oscillator 
strength,^  the  geometry  of  nanocrystals  and  its  effect  on  the  electronic  structure,  2c,3,4  ^nd  the 
location  of  electron-hole  pairs  in  real  space.^’^ 

Recent  photoluminescence  experiments  suggest  strong  resonance  mixing  between  surface 
and  core  states.^^  Also,  it  is  well  known  that  surface  states  arise  from  loss  of  periodicity  at 
crystal  surfaces.'^*’’  These  observations  raise  questions  concerning  electron  and  hole 
locations,  possible  preferred  sites  in  the  excited  state,  energetic  position  of  surface  states  (caused 
by  dangling  bonds  /  incomplete  capping  on  the  surface)  and  their  possible  mixing  with  other 
states.  If  an  electron-hole  pair  is  separated  spatially  in  the  excited  state,  then  a  change  in  electron 
population  from  ground  to  excited  state  will  be  observed,  i.e.  a  significant  dipole  moment 
change.  Therefore,  by  creating  an  incomplete  capping  site  on  the  cluster  surface  and  analyzing 
the  change  in  electron  population,  we  can  further  understand  these  issues. 

Theoretical  models  such  as  the  effective  mass  model,^^’^'^  the  tight-binding  model,  time- 
dependent  tight-binding  approach, and  the  empirical  pseudopotential  methodic  have  been 
employed  to  understand  the  electronic  properties  of  nanocrystals. 

The  effective  mass  model  seems  adequate  for  larger  particles  but  overestimates  the  band  gap 
for  smaller  ones;  it  considers  a  nanocrystal  as  a  particle  in  a  spherically  symmetric  well  with  an 
infinite  wall  on  the  particle  surface.  Recently,  one  group  modified  the  effective  mass  model  by 
taking  account  of  the  valence  band  degeneracy  and  the  conduction  band  nonparabolicity.^^l 
Results  are  in  good  agreement  with  experiment  for  rather  large  particles.^^^  However,  the  surface 
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localization  of  electron  and/or  hole  has  been  ignored  by  imposing  a  boundary  condition  such  that 
both  the  electron  and  hole  wavefunctions  vanish  at  the  particle  surface. 

The  empirical  pseudopotential  methodic  yields  better  agreement  with  experiment.  This 
model  describes  the  nanoparticle  as  a  bulk  material,  allowing  the  use  of  the  periodic  boundary 
conditions  (or,  introducing  the  reciprocal  lattice  vector  G).  The  particle  interface  is  taken  as  an 
^nite  wall.  Recent  work  modifies  the  form  factors  of  both  zinc -blende  and  hexagonal 
structures  to  reproduce’ the  bulk  crystal  field  splitting  at  the  top  of  the  valence  band.  13  However, 
the  problem  of  surface  effects  is  still  not  addressed. 

The  band  gap  energies  calculated  by  the  tight-binding  scheme^O  and  the  time-dependent 
tight-binding  approach  are  also  in  reasonable  agreement  with  experiment  for  intermediate  to 
larger  particles.  For  the  latter  approach,  ^  ^  surface  effects  have  been  considered  by  covering  the 
surface  with  dangling  bonds  and  reparameterizing  the  interactions  between  surface  atoms  and 
ligand  atoms.  The  electron-electron  repulsion  is  included  within  the  electron  hole  interaction. 
However,  the  Huckel-type  Hamiltonian,  parameterization  from  the  bulk  band  structure,  and 
consideration  of  only  nearest  neighbor  interactions  make  it  difficult  for  these  models  to  predict 
the  band  gap  over  a  broad  size  range. 

The  semiempirical  methods  of  electronic  structure  calculation  have  been  widely  employed  by 
researchers  in  many  areas  to  understand  the  electronic  properties  of  molecules.  Such 

calculations  tire  able  to  describe  efficiently,  intuitively,  and  quite  accurately  the  electronic 
properties  of  ground  and  excited  states.  Our  goal  here  is  to  understand  the  nature  of  the 
electronic  transitions  for  CdS  nanoclusters,  the  generation  of  an  optically  excited  electron  in  a 
cluster,  the  effects  of  a  defect  site  on  the  cluster  surface,  and  the  effects  of  solvation  on 
absorption  spectra  by  using  the  complete  neglect  of  differential  overlap  (CNDO)  model 
hamiltonian  and  the  configuration  interaction  singles  (CIS)  approximation.  For  comparison  with 
CNDO  CIS  results,  the  extended  Huckel  method  has  also  been  used  to  calculate  the  orbital 
energies.  Existing  crystal  data  are  used  to  determine  the  geometries.  Due  to  limitations  of 
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computer  resources  and  parameter  availability,  only  CdS  molecules/clusters  of  three  sizes  are 
examined  here. 

In  section  II,  we  summarize  the  models  and  analyze  parameters  used.  Section  III  discusses, 
in  order,  the  orbital  energies  from  both  EHT  and  CNDO  calculations,  the  simulated  optical 
spectra  and  assignment  of  the  lowest  energy  band,  the  effects  of  a  defect  site  located  at  one  Cd 
atom,  and  the  effects  of  solvation.  Section  IV  gives  some  conclusions. 

11.  Computational  Details 
A.  Hamiltonians  and  states 

Two  semiempirical  models  are  employed  here:  the  one-electron  extended  Huckel  model 
(EHT)  and  the  complete  neglect  of  differential  overlap  (CNDO)  model,  that  includes  electronic 
repulsions.  The  resonance  integral,  Hij,  in  the  EHT  calculation  is  evaluated  via  the  Wolfsberg- 
Helmholz  equation  with  K  equal  to  1.75.  To  compute  experimental  spectra  more  accurately, 
the  CNDO  model  with  monoexcited  Cl  or  CIS,  has  been  employed.  The  computational  tool  is 
the  CNDO  routine  in  the  ZDIDO  program  of  Zemer  and  co-workers. The  matrix  elements  of 
the  Fock  operator  are: 

Fnn  =  U^l^-t■(PAA-l/2Pp^l)TAA+  I  (PbbT^A-Vab)  (D 

B(?iA) 

Fjxv  =  6^  Sjiv  ■  l/2P|iv  Tab  (1  ^  V  (2) 

The  atomic  matrix  elements  of  the  one-electron  Hamiltonian  U|x|j,  are  determined  from  the 

corresponding  ionization  potential.  The  one-center-two-electron  Coulomb  integral  Yaa  has 
three  different  values:  ,  and  for  transition  metal  Cd^O  and  only  one  value 

for  others.  For  the  optical  spectra,  the  Mataga-Nishimoto  formula  is  employed  to  evaluate  the 
two-center-two-electron  integrals  Yab  (including  7^  ,  7^  ,  and  7^  ):2^ 
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(3) 


a  +  Rab 

7 

a  = - = -  (4) 

'-'aa  +  '^bb 

The  bonding  parameter  in  eqn(2)  is  calculated  also  by  the  Wolfsberg-Helmholz 

approximation  with  appropriate  interaction  factors  parameterized  for  a  and  n  bonds. -2  Vab 
S|j,v  are  calculated  as  described  in  ref.  16(b). 

The  lowest  transition  energies  and  oscillator  strengths  are  calculated  using  CIS,  including 
less  than  700  singly  excited  configurations  of  valence  orbitals.  To  assure  that  the  number  of 
singly  excited  configurations  is  large  enough  to  converge,  tests  have  been  done  on  the  calculated 
transition  energies  as  a  function  of  the  number  of  configurations  at  21 1.3  nm,  201.8  nm,  and 
223.6  nm  for  [Cd(SPh)4]2-  ([Cdl]),  [Cd2(SPh)6]2-  ([Cd2]),  and  [Cd4(SPh)io]2-  ([Cd4]) 
respectively.  Results  show  that  enough  configurations  are  included:  transition  energies  are 
unchanged  to  within  100  crn'l  or  so  (Appendix  A).  The  state  dipole  moments  are  evaluated  as 
expectation  values  of  the  dipole  moment  operator.  For  the  elements  C,  H,  S,  the  default 
parameters  are  employed  without  further  refinement  to  fit  the  experimental  data  (Appendix  B). 
For  Cd  the  Slater-Condon  F®  (Taa)  integral  is  taken  from  ref.  23.  The  geometries  of  [Cdl], 
[Cd2],  and  [Cd4]  used  in  both  EHT  and  CNDO  calculations  are  from  crystal  structure 
determinations  (Scheme  l).24-26  Experimental  optical  absorption  spectra  are  from  the 

literature.25.27 

B.  Lineshapes 

CIS  calculation  results  give  discrete  transition  energies  and  oscillator  strengths  while  an 
experimental  absorption  spectrum  yields  a  curve  of  the  wavelength-dependent  extinction 
coefficient  function.  To  compare  the  calculated  spectrum  with  empirical  bands  and  also  to  take 
account  of  the  Franck-Condon  broadening,  it  is  convenient  to  use  the  Lorentzian  lineshape  to 
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simulate  each  calculated  electronic  transition.-^  The  normalized  Lorentzian  lineshape  with 
extinction  coefficient  function,  £’(v),  centered  at  Vq  can  be  expressed  as: 


e’fv)  = 


^max 

K 


(r/2) 

(V-Vo)2  +  (r/2)2 


(5) 


where  Smax  is  the  extinction  coefficient  at  the  band  maximum,  Vq,  and  T  is  the  band  width  at  half 
maximum  (FWHM).  The  oscillator  strength  for  this  band  is 


f(Vo)  =  k  j  £’(v)  dv 


(6) 


where  k  is  a  constant  equal28.29  jq  4.33  x  10"^  M’cm^. 

The  integral  on  the  right  hand  side  of  eqn(6),  called  the  integrated  extinction  coefficient, 
extends  over  the  entire  band.  To  correlate  f  with  emax.  the  Lorentzian  lineshape  function 
(eqn(5))  can  be  fed  into  eqn(6): 


f(Vo)  =  k  je’(v)dV  =k£niax 


or 


t'max  — 


Replacement  of  £max  in  eqn(5)  with  eqn(8)  yields 


Vo  k  (v-Vo)2  +  (r/2)2 


=  1.155  X  108 


f(vo)  r 


(V-Vo)2  +  (r/2)2 


(7) 

(8) 
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£’  (V )  ,  in  units  of  M'  *  cm" ' ,  denotes  the  molar  extinction  coefficient^^  associated  with  the 

transition  occurring  at  wavenumber  Vq.  Since  there  are  transitions  at  various  energies,  the 
calculated  spectra  are  obtained  by  superposition  of  each  £’^  fv)  band: 


e(v)  =  I 

Vo 


=  X 

Vo 


1.155  X  108  - - 

(V-Vo)2  +  (r/2)2 


(10) 


e(v)  is  the  simulated  spectrum  used  to  compare  with  the  experimental  one.  T  is  adjusted  to  fit 
experimental  data. 

C.  Electron  populations 

Electron  density  plots  are  generally  pictorially  inadequate  to  show  how  electron  density 
changes  during  an  excitation  process  for  unsymmetrical  molecules  such  as  [Cdl],  [Cd2],  and 
[Cd4].  Therefore,  estimation  of  the  changes  in  electron  population  for  each  element  from 
ground  to  excited  state  is  needed  to  understand  the  nature  of  each  transition.  We  can  introduce  a 
very  approximate  way  to  evaluate  the  changes  in  population.  Assume  that  an  electron  in  the  pth 
molecular  orbital  of  a  compound  in  the  ground  state  is  excited  to  the  vth  virtual  orbital.  Then  the 
total  change  in  the  total  number  of  electrons,  assuming  zero  differential  overlap,  is 


0  =  An  =  nex  -  ng 


A.O.OCC 

=  I  I  I 

A  i  o 


laA 


ex 

SaA 


A.Oocc 

III 

A  J  y 


o* 

C® 

J  ya 


c® 

j  yA 


(11) 


Here  A,  c,  and  i  label  atoms,  molecular  orbitals  and  atomic  orbitals  respectively.  For  a  single 
excitation,  common  occupations  cancel,  and  therefore 
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(12) 


0  =  An  =  X 
A 


A.O. 

I(C 


I 


ex*  ex 
i  vaS'  vA 


a*  o 

/  |IA  /  |IA^ 


All  terms  cancel  except  those  M.O.  involved  in  the  electronic  transition  |j,  ->  v.  Since  each 
excited  state  consists  of  several  singly  excited  configurations,  the  resulting  total  change  in 
Mulliken  population  is  obtained  by  appropriately  weighted  contributions  from  all  of  the 
participating  configurations. 


0  =  AN  =  -  NS 


config. 

=  I  I 

A  k(|j.-»v) 


A.O. 

yic®’"  c®"'  -c^  c®  i 

4  ''  vAkS'  vAk  S'  pAk^j  ^cAk'' 


(13) 


* 

dj^  dj^  is  the  probability  for  configuration  k  to  occur  and  is  determined  in  CIS  calculations. 

Here  only  those  configurations  with  Idkl  greater  than  0. 1  are  considered.  Therefore  for  element 
A,  the  change  in  Mulliken  population  from  the  ground  state  to  an  excited  state  can  be  expressed 
as: 


config.  A.O. 

ANa=  I  I 

k(|i->v)  I 


ex*  ex 
SvAkSvAk 


S°^Ak^r|lAk^ 


(14) 


ANa  has  the  value  between  0  and  1.  Thus  ANa  x  100%  represents  the  percentage  of  one 
electron  that  is  gained  by  atom  A  during  an  excitation  process. 

D.  Solvent  effects 

The  effect  of  solvent  is  examined  by  applying  the  self-consistent  reaction  field  (SCRF) 
method.30  To  exclude  the  volume  of  counter  ions,  the  size  of  the  solute  molecules  is  evaluated 
through  the  solute  volume  density  obtained  from  crystal  data. 

E.  Defects 
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The  effects  on  electronic  spectra  due  to  incomplete  capping  of  the  cluster  surface  are 
explored  by  removing  one  terminal  ligand  from  the  studied  molecule.  All  defect  molecules 
retain  the  same  geometries  and  no  surface  reconstruction  is  considered. 


III.  Results  and  Discussion 
A.  Frontier  orbitals  of  CdS  molecules 

For  bulk  CdS  semiconductor  materials,  the  frontier  orbitals  consist  of  S  atoms  as  the  valence 
band  and  Cd  atoms  as  the  conduction  band.  On  the  other  hand,  for  molecular-sized  CdS 
compounds,  the  loss  of  periodicity  or  variations  in  the  coordinating  environment  causes  different 
redox  properties  for  the  same  kind  of  atoms  in  a  molecule.  That  results  in  removal  of 
degeneracy  occurring  in  both  Cd  and  S  bands,  and  therefore  different  bands  of  optical  transitions 
arise.  For  example,  [Cd2]  and  [Cd4]  clusters  contain  both  bridging  and  capping  thiolate 
ligands:  The  S  atom  of  a  bridging  thiolate  has  two  neighboring  Cd  atoms,  but  the  S  atom  of  a 
terminal  thiolate  shares  its  electrons  with  only  one  Cd  atom.  Thus  the  terminal  sulfur  is 
anticipated  to  be  more  electron-donating  than  the  bridging  one.  Both  EHT  and  CNDO 
calculations  show  that,  even  though  the  terminal  and  bridging  sulfurs  have  distinct  redox 
properties  and  the  terminal  sulfurs  do  lie  higher  in  energy  than  bridging  ones  in  both 
calculations,  the  energy  difference  is  small.  Figures  1  and  2  show  the  orbital  energies  from  EHT 
and  CNDO  respectively.  The  two  models  give  different  results.  In  Figure  1,  EHT  results  show 
that  orbitals  sitting  between  HOMO  and  the  empty  orbitals  centered  at  Cd  5s  orbitals  (with 
energies  greater  than  -2  eV)  are.primarily  comprised  of  antibonding  K  (k*)  orbitals  localized  on 
thiolates.  In  other  words,  Cd  5s  orbitals  nearly  do  not  mix  with  thiolate  k*  orbitals.  Thus,  from 
the  EHT  picture,  the  lowest  energy  transitions  will  be  only  from  sulfur  to  k*  orbitals  within 
thiolates;  thiolate  k  to  Cd  5s  transitions  will  require  much  higher  energies  than  the  energies  for 
intraligand  charge  transfer  transitions.  If  we  generate  an  EHT  spectrum,  the  stronger  intraligand 
transitions  will  lie  lower  in  energy  than  the  weaker  thiolate  to  Cd  transitions.  This  result  is 
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obviously  not  reasonable:  the  band  shape  features  of  experimental  spectra  show  that  the  stronger 
band  has  higher  energy  than  the  weaker  one. 

On  the  other  hand,  from  the  results  of  CNDO  calculations  in  Figure  2,  the  n*  orbitals  on 
thiolates  do  mix  with  Cd  5s  orbitals.  As  usual  with  SCF  calculations,  an  artificially  large  gap 
between  HOMO  and  LUMO  appears.  Also  note  that  in  Figure  2,  the  Cd  5s  virtual  orbital 
energies  stay  about  the  same  from  [Cdl]  to  [Cd4]  but  the  energies  of  the  highest  occupied 
molecular  orbitals  decrease  as  the  size  of  the  molecule  increases.  Evidently  the  electron-electron 
repulsion  considered  in  the  CNDO  calculation  pushes  up  the  empty  thiolate  7t*-based  molecular 
orbitals  to  about  the  same  energy  level  as  the  Cd  5s  orbitals.  Although  for  larger  CdS  clusters, 
the  ratio  of  the  number  of  stabilizing  ligands  (e.g.,  thiolates)  versus  the  number  of  sulfur  atoms  is 
small,  a  recent  study  has  shown  that  the  ligands  are  essential  to  stabilize  smaller  clusters.^ 
Therefore,  the  role  of  ligands  cannot  be  ignored  for  small  clusters. 

B.  Electronic  Absorption  Spectra  and  Changes  in  Dipole  Moment. 

Since  they  account  for  incorrect  virtual  orbital  energies,  the  CNDO  CIS  calculations  yield 
lower  excitation  energies  in  comparison  with  orbital  energy  differences  only.  The  nature  of  the 
lowest  energy  transitions  can  be  understood  by  performing  an  analysis  of  the  changes  in 
Mulliken  population  associated  with  various  transitions  (Figures  3-5  and  Tables  2-4). 

The  CNDO  Cl  results  show  that  the  calculated  spectra  well  reproduce  the  main  features  of 
the  experimental  bands.  The  experimental  UV-VIS  spectrum  for  [Cdl]  peaks  at  282  nm  with  a 
broad  band  extending  from  300. to  below  200  nm  (Table  1).27  This  band  was  originally  assigned 
to  ligand  to  metal  charge  transfer  (LMCT)  transitions.^^  In  the  CNDO  Cl  calculated  spectrum, 
the  lowest  energy  band  appears  in  the  correct  wavelength  range  (3 10  to  215  nm),  but  is 
dominated  on  the  red  side  by  intraligand  transitions  (see  the  lower  plot  in  Figure  3).  The 
calculated  strong  LMCT  transitions  occur  on  the  blue  side  of  the  spectrum  (Table  2). 

The  optical  absorption  spectrum  for  [Cd2]  has  a  maximum  at  270  nm  and  has  been 
previously  attributed  to  both  intraligand  and  thiolate  to  Cd  transitions  (Table  1).25  In  Figure  4, 
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the  stronger  transitions  within  the  band  centered  at  240  nm  clearly  include  both  LMCT  and 
intraligand  transitions,  but  with  a  greater  contribution  from  the  former.  Toward  the  red  side  of 
the  band,  terminal  thiolates,  and  to  a  lesser  extent,  bridging  thiolates,  contribute  charge  transfer 
(CT)  intensity.  The  differences  in  degree  of  participation  are  clearly  understandable  in  light  of 
the  slightly  higher  orbital  energies  for  the  terminal  versus  bridging  ligands. 

The  [Cd4]  experimental  spectrum  exhibits  a  maximum  at  249  and  an  overlapping  band  at 
275  nm,  with  tails  extending  from  320  to  below  200  nmP  The  combined  band  has  previously 
been  assigned  as  a  ligand  to  metal  charge  transfer  transition.^^  As  shown  by  Figure  5,  the 
CNDO  Cl  calculations  yield  transitions  at  higher  energies  than  found  experimentally.  However, 
the  overlapping  two-band  feature  is  observed  computationally.  From  Table  4,  intraligand 
excitations  clearly  must  play  the  major  role  for  both  absorption  bands,  since  the  LMCT 
transitions  have  very  small  oscillator  strengths.  Therefore,  for  [Cdl],  [Cd2],  and  [Cd4],  the 
thiolate  n*  and  Cd  5s  orbitals  mix  to  make  up  the  low-lying  excited  states  and  the  terminal  and 
bridging  thiolates  contribute  mostly  to  the  ground  state.  The  increase  in  extinction  coefficient 
with  increasing  molecular  size  is  due  to  increasing  number  of  atoms  in  a  larger  cluster  which 
provides  more  available  states  for  transitions,  rather  than  greater  oscillator  strength  for  each 
individual  transition.  The  intraligand  transitions  dominate  strongly  in  the  [Cd4]  species. 

An  “anomalous”  blue  shift  with  increasing  cluster  size  is  observed  both  experimentally  and 
theoretically  (Figures  3  -  5).  This  shift  is  consistent,  however,  with  the  systematic  decrease  in 
HOMO  energy  with  increasing  cluster  size  (Figure  2).  The  decreasing  HOMO  energy  effect  can 
be  explained  by  the  Jahn-Teller.effect,3l  which  says  that  the  ground  state  of  a  molecule  is 
stabilized  by  undergoing  geometry  perturbations  in  order  to  remove  degeneracies  in  the 
unperturbed  ground  state  and  thus  lower  the  total  energy  of  the  molecule.32  Consistent  with  this 
idea,  the  [Cdl]  molecule  lacks  a  Tj  point  group  symmetry.  The  bond  lengths  between  Cd  and  S 
atoms  vary  from  2.517  to  2.546  Similarly  the  [Cd2]  molecule  has  short  terminal  bond 
lengths  and  long  bridging  bond  distances  between  Cd  and  S  atoms.  [Cd4]  also  experiences 
geometry  distortions:  four  different  terminal  Cd-S  bond  lengths  (ranging  from  2.460  to  2.480 
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and  six  different  bridging  Cd-S  bond  lengths  (ranging  from  2.533  to  2.603  A)  are  observed 
in  the  [Cd4]  molecule.-^  Thus  [Cd4]  has  no  symmetry’  elements  (  C)  point  group).  These 
geometry  distortions  in  [Cdl],  [Cd2],  and  [Cd4]  systematically  lower  the  energy  levels  of 
thiolates  and  thus  blue  shifts  in  the  optical  spectra  are  seen. 

In  contrast  to  [Cdl],  [Cd2],  and  [Cd4],  [S4Cdio(SPh)i6]^'  ([Cdl0])36  does  possess  several 
sets  of  Cd  and  S  atoms  that  are  crystallographically  equivalent  to  other  atoms  in  the  same 
molecule.  These  equivalent  sets,  originating  from  the  crystal  symmetry,  bring  together 
discrete  levels  to  yield  bands;  the  bandwidths,  which  are  determined  by  the  number  of 
participating  atoms  or  moieties,  play  a  major  role  in  determining  the  energies  of  transitions  from 
the  "valence  band"  to  the  "conduction  band". 

For  bulk  materials,  surface  relaxation  and  surface  reconstruction  are  the  two  dominant  factors 
defining  whether  surface  states  are  significant.^*  As  the  size  of  a  material  decreases,  the  ratio  of 
surface  atoms  to  "core"  atoms  obviously  increases.  In  the  limit  of  very  small  particle  sizes  (i.e. 
clusters),  therefore,  surface  states  and  geometry  distortions  may  well  dominate  the  electronic 
structure. 

The  question  of  how  electron  density  is  redistributed  upon  optical  excitation  clearly  is  related 
to  the  issue  of  how  the  dipole  moment  is  changed  upon  conversion  of  the  clusters  from  ground 
state  to  an  excited  state.  For  example,  in  Table  2,  transitions  for  [Cdl]  at  294.5  nm  and  266.9 
nm  can  both  be  characterized  as  intraligand  charge  transfer.  However,  their  dipole  moment 
changes  differ  (5.3  D  vs.  14.1  D).  Further  study  on  the  population  changes  shows  that  the  294.5 
nm  transition  mainly  occurs  on  one  ligand  while  the  266.9  nm  transition  involves  charge  transfer 
from  one  ligand  to  another  -  thus  the  larger  dipole  moment  change  for  the  latter.  On  the  other 
hand,  the  small  net  dipole  moment  change  at  227  nm  (Table  2)  reflects  simultaneous  charge 
transfer  from  sulfur  to  both  cadmium  and  benzene.  This  results  in  partial  cancellation  of  changes 
in  dipole  moment.  Small  dipole  moment  changes  are  also  encountered  for  transitions  involving 
[Cd2]  (Table  3).  Here,  however,  an  additional  factor  may  come  into  play.  The  cadmium  atoms 
in  [Cdl]  and  [Cd4]  molecules  do  not  have  a  particular  point  group  symmetry  associated  with 
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them.  The  [Cd2]  molecule,  on  the  other  hand,  possesses  a  reflection  center  (i).  This  center 
inhibits  electron  localization.  Since  the  hole  is  delocalized,  localization  of  an  excited  electron 
causes  a  significant  dipole  moment  change,  while  delocalization  results  in  a  small  change  in 
dipole  moment. 

C.  Effects  of  Defect  Sites  on  Optical  Spectra. 

Experimental  studies  on  surface  states  and  exciton  traps  of  quantum-confined  particles 
suggest  that  surface  electronic  levels  mix  with  interior  levels.^^-^*’  Surface  states  can  arise  either 
from  differences  of  redox  properties  of  surface  moieties  from  the  core,  which  implies 
delocalization  of  surface  states,  or,  from  incomplete  coverage  of  the  particle  surface,  which  may 
imply  localization  of  surface  states.39-45  jjgre  we  try  to  mimic  the  effects  of  incomplete  surface 
capping  by  removing  one  terminal  thiolate  from  [Cd2]  and  [Cd4].  CNDO  results  show  that  a 
dangling  bond  on  Cd  atom  creates  a  distinct  low-lying  Cd  5s  empty  orbital  which  acts  as  an 
electron  acceptor/trap  upon  an  electronic  excitation.  In  Figures  6  and  7,  new  excitations  from 
both  terminal  and  bridging  thiolates  to  the  defect  cadmium  generate  new  transitions  with  small 
oscillator  strengths  at  the  red  end  of  a  spectrum  (cf.  Figures  4  and  5).  Although  terminal 
thiolates  play  a  major  role  as  electron  donors,  there  again  are  no  significant  energy  differences 
between  terminal  and  bridging  thiolates.  In  addition,  as  shown  in  the  lower  plots  in  both  figures, 
these  transitions  have  a  strong  ligand  to  metal  charge  transfer  character.  As  one  might  expect,  a 
significant  dipole  moment  change  associated  with  each  "LMCT"  excitation  is  obvious  (Tables  5, 
6). 

Computational  limitations  preclude  similar  calculations  on  larger  clusters  such  as 
[S4Cdio(SPh)i6l2',^^  [Si3Cd2o(SPh)22]^','^^  and  Si4Cd32(SPh)36;‘^^  however,  the  long  red  tails 
and  the  corresponding  significant  dipole  moment  changes  found  for  the  defect  structures  (but 
absent  in  Figures  2-4),  suggest  that  if  exciton  trapping  is  due  to  incomplete  capping  on  the 
surface,  then  exciton  localization  and  a  large  dipole  moment  change  (perhaps  observable  by 
Stark  spectroscopy)  should  generally  be  expected.  We  note  that  the  degree  of  dipolar  character 
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of  excited  states  for  ver>'  large  clusters  is  still  controversial  experimentally. 5-6  On  the  other 
hand,  the  small  oscillator  strengths  associated  here  with  large  A|i  transitions  clearly  suggest  that 
discrete  strong  optical  transitions  from  the  delocalized  ground  state  to  such  trapped  states  can  not 
occur,  at  least  if  the  surface  trap  states  are  ascribed  to  incomplete  coordination  of  the  surface. 

Due  to  the  poor  overlap  between  the  ground  state  and  the  Cd-localized  trap  state,  it  is  unlikely 
that  light  can  be  emitted  efficiently  (i.e.,  both  rapidly  and  with  high  quanmm  efficiency)  upon 
relaxation  from  the  excited  Cd  trap  state.  This  is  supported  by  recent  experimental  studies.-^^ 
Nevertheless,  basic  questions  of  the  nature  of  surface  states  and  correlation  of  an  exciton  trap 
with  the  surface  redox  properties,  incomplete  coverage  on  the  surface,  and  even  the  packing 
faults  of  clusters  are  still  uncertain  and  require  more  study. 

D.  Solvent  Effect 

By  using  the  SCRF  model  to  examine  the  effect  of  solvation  of  clusters  by  acetonitrile,  we 
see  shifts  in  optical  transition  energy  ranging  from  several  cm’l  to  (at  most)  ca.  400  cm'l. 
Therefore  the  simulated  spectra  suffer  no  significant  change  upon  inclusion  of  solvation. 
According  to  the  SCRF  model,  the  shift  in  energy  due  to  solvation  will  be  proportional  to  the 
square  of  dipole  moment  and  inversely  proportional  to  the  cube  of  solute  radius,  Since 

even  small  clusters  have  rather  large  radii  by  molecular  standard  (5.4,  6.3,  and  7.4  A  respectively 
for  [Cdl],  [Cd2],  and  [Cd4] ),  even  a  change  in  dipole  moment  as  great  as  lOD  will  decrease  the 
transition  energy  by  only  a  few  hundred  wavenumbers.  The  same  situation  is  observed  in 
calculations  on  defect  molecules.  The  SCRF  model  thus  predicts  only  small  absorption  energy 
shifts.  However,  it  is  not  clear  that  the  SCRF  model  is  adequate  for  these  clusters.  Of  particular 
concern,  it  is  the  assumption  of  a  strictly  spherical  solvent  cavity  with  no  specific  local 
interaction  between  solute  and  solvent.  It  neglects  specific  effects  such  as  those  brought  about 
by  solvent  coordination  at  trap  state.  Experimental  studies  show  that  coordinating  solvents  can 
stabilize  nanocrystallites  and  increase  solubility;^^  both  observations  suggest  roles  for  local 
interactions  between  particle  and  solvent.  Furthermore  small  CdS  clusters  have  a  decidedly  non- 
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spherical  shape.  36, 46,47,49  this  early  stage,  even  an  experimental  consensus  on  solvent  effects 

is  still  lacking.3,46,50 


IV.  Conclusions 

The  simulated  spectra  in  the  CNDO  calculations  well  reproduce  the  main  features  of 
experiments.  The  lowest  energy  band  is  made  up  of  transitions  with  both  LMCT  (thiolates  to  Cd 
5s)  and  intraligand  charge  transfer  character.  The  blue  shift  with  increasing  cluster  size, 
observed  both  experimentally  and  theoretically,  can  be  explained  by  the  Jahn-Teller  effect  which 
emphasizes  the  importance  of  geometric  distortions  to  the  electronic  structure  of  cluster.  Also 
accompanying  these  distortions  are  reduction  in  symmetry,  resulting  in  dimension  of  orbital 
degeneracy  and  enhanced  localization.  If  the  electron  is  localized  on  one  site  (either  on  a  single 
Cd  atom  or  a  thiolate)  in  an  excited  state,  a  significant  change  in  dipole  moment  from  the  ground 
to  the  excited  state  results.  However,  if  a  molecule  possesses  an  i  symmetry  center,  vanishing 
changes  in  dipole  moment  are  expected. 

Removal  of  one  terminal  ligand  creates  a  defect  site  on  a  Cd  atom.  This,  in  turn,  yields 
several  weak  transitions  associated  with  excitations  from  thiolate  ligands  to  the  incompletely 
capped  Cd.  These  transitions  collectively  appear  as  a  long  tail  on  the  red  side  of  the  spectrum. 
Significant  changes  in  dipole  moment,  attributed  to  localization  of  an  excited  electron,  are 
associated  with  these  transitions.  However,  the  small  oscillator  strengths  associated  with  these 
transitions  indicate  poor  overlap  between  the  ground  and  the  trap  states,  which  suggests  that  it  is 
unlikely  that  the  Cd  atoms  with  incomplete  capping  can  serve  as  both  as  surface  trap  states  and 
efficient  light-emitting  states. 

Solvent  effects  have  also  been  examined  by  using  the  SCRF  model.  Due  to  the  large 
molecular  radii,  no  significant  changes  in  simulated  spectra  are  observed. 
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Appendix  A 

To  assure  that  the  Cl  calculation  is  converged  as  a  function  of  the  number  of  excitations, 
plots  of  transition  energy  versus  number  of  configurations  are  made  for  three  transitions  chosen 
from  the  blue  side  of  each  [Cdl],  [Cd2],  and  [Cd4]  spectrum. 

From  Figures  A1  and  A2,  the  number  of  configurations  included  in  Cl  calculations  (691  and 
677  for  [Cdl]  and  [Cd2]  respectively)  place  the  calculations  in  the  flat  regions  of  the  curves, 
indicating  that  convergence  has  occurred. 

In  Figure  A3,  although  the  number  of  atomic  orbitals  increases  rapidly  with  increasing 
cluster  size,  the  transition  energy  itself  becomes  stable  if  the  number  of  configurations  is  greater 
than  650. 

Appendix  B 

The  default  ionization  potentials  (U^x),  bonding  parameters  (Pa),  and  one-center-two- 
electron  Coulomb  integrals  (Taa)  for  H,  C,  and  S  are  taken  as^^: 


Uss  (eV) 

UDo(eV) 

PA(eV) 

v  * 
»AA 

H 

-13.06 

-3.23 

-12 

12.85 

C 

-19.65 

-11.13 

-17 

11.11 

S 

-22.76 

-12.08 

-35 

10.09 

*  Yss  -  Ypp  —  Ysp 
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Figure  captions 

Figure  1.  Orbital  energies  in  EHT  calculations  for  [Cd(SPh)4]2-  (top),  [Cd2(SPh)6p‘  (middle), 
and  [Cd4(SPh)ioF'  (bottom).  The  solid  lines  represent  occupied  molecular  orbitals.  The  dash 
lines  are  from  virtual  orbitals. 

Figure  2.  Orbital  energies  in  CNDO  calculations  for  [Cd(SPh)4]2-  (top),  [Cd2(SPh)6p'  (middle), 
and  [Cd4(SPh)io]^‘  (bottom).  Layout  as  in  Figure  1. 

Figure  3.  (upper)  CNDO  oscillator  strengths  (virtual  lines)  and  simulated  spectrum  (solid  curve) 
for  [Cd(SPh)4]2-  ([Cdl]);  inset,  experimental  absorption  spectrum  from  ref.  23.  (lower)  Changes 
in  electron  population  from  the  ground  to  various  excited  states.  Only  dominant  transitions  are 
selected.  Solid  line/full  circle,  Cd;  dash  line/square,  S;  dotted  line/triangle,  C. 

Figure  4.  CNDO  oscillator  strengths,  simulated  spectrum,  and  changes  in  electron  population  for 
[Cd2(SPh)6]2‘  ([Cd2]).  Layout  as  in  Figure  3. 

Figure  5.  CNDO  oscillator  strengths,  simulated  spectrum,  and  changes  in  electron  population  for 
[Cd4(SPh)ioF'  ([Cd4]).  Layout  as  in  Figure  3. 

Figure  6.  CNDO  oscillator  strengths,  simulated  spectrum,  and  changes  in  electron  population  for 
[Cd2(SPh)5]'.  Layout  also  as  in  Figure  3.  Note  that  more  transitions  with  small  oscillator 
strength  appear  at  the  red  side  of  the  spectrum  and  they  have  a  character  of  S  and  n  to  Cd  charge 
transfer.  See  text  also. 

Figure  7.  CNDO  oscillator  strengths,  simulated  spectrum,  and  changes  in  electron  population  for 
[Cd4(SPh)9]*.  Layout  as  in  Figure  3.  Note  that  transitions  with  small  oscillator  strength  at  the 
red  side  of  the  spectrum  correspond  to  the  charge  transfer  from  S  and  it  to  Cd.  See  text  also. 
Figure  Al.  Number  of  singly-excited  configurations  versus  transition  energy  for  [Cdl].  As 
more  than  300  configurations  are  included  in  the  Cl  calculation,  the  transition  energy  curve 
becomes  a  flat  line. 

Figure  A2.  Number  of  singly-excited  configurations  versus  transition  energy  for  [Cd2].  The  flat 
region  starts  in  which  the  number  of  configurations  is  greater  than  350. 

Figure  A3.  Number  of  singly-excited  configurations  versus  transition  energy  for  [Cd4]. 
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Table  2.  Transilion  energies",  oscillator  strengths,  changes  in  Mulliken  population^,  and  dipole  moment  changes^-’  of 
[Cd(SPh)4]2-. 


b.  It  is  evaluated  by  weighted  Mulliken  population  differences  of  molecular  orbitals  involved  in  an  electronic 
transition.  Also  see  text. 

c.  It  is  the  dipole  moment  change  from  the  ground  state  to  a  singly  excited  state. 


Table  3.  Transition  energies^,  oscillator  strengths,  changes  in  Mulliken  population*’,  and  dipole  moment  changes‘^  of 
[Cd2(SPh)6]2-. 


Table  4.  Transition  energies®,  oscillator  strengths,  changes  in  Mulliken  population*’,  and  dipole  moment  changes^^  of 
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Superscripts  a,  b,  and  c  as  in  Table  2. 

(p,  d,  b,  t )  denotes  (Cd  with  complete  coordination,  Cd  with  incomplete  coordination,  bridging  ligand,  terminal  ligand). 


Table  6.  Transition  energies^,  oscillator  strengths,  changes  in  Mulliken  population*’,  and  dipole  moment  changes*^  of 
[Cd4(SPh)9]-. 
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